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FRACTIONAL DIFFUSION EQUATIONS AND PROCESSES WITH 
RANDOMLY VARYING TIME^ 

By Enzo Orsingher and Luisa Beghin 

"Sapienza" Universitd di Roma 

In this paper the solutions Ui, — Uu{x,t) to fractional diffusion 
equations of order Q <v <2 are analyzed and interpreted as densities 
of the composition of various types of stochastic processes. 

For the fractional equations of order v = r^, n > 1, we show that 
the solutions Mi/2" correspond to the distribution of the n-times it- 
erated Brownian motion. For these processes the distributions of the 
maximum and of the sojourn time are explicitly given. The case of 
fractional equations of order j/ = n > 1, is also investigated and 
related to Brownian motion and processes with densities expressed 
in terms of Airy functions. 

In the general case we show that Uv coincides with the distribution 
of Brownian motion with random time or of different processes with 
a Brownian time. The interplay between the solutions and stable 
distributions is also explored. Interesting cases involving the bilateral 
exponential distribution are obtained in the limit. 

1. Introduction. Time- fractional equations of the form 



u 



(1-1) w = >-"ii^^- 

for < < 2, have been studied by a number of authors since the 1980s: see, 
for example, Wyss (1986), NigmatuUin (1986), Schneider and Wyss (1989), 
Mainardi (1995a, 1996) and, more recently, NigmatuUin (2006), Angulo et 
al. (2000, 2005). Hyperbolic fractional equations similar to (1.1) have been 
analyzed, for example, by Engler (1997). 

For exhaustive reviews on this topic, also consult Samko, Kilbas and 
Marichev (1993) and Podlubny (1999). 
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For interesting applications of fractional equations to physical problems 
see, for example, Saichev and Zaslavsky (1997), Nigmatullin et al. (2007), 
Angulo et al. (2005). 

Fractional diffusion equations of order 1 < < 2 emerge in the study of 
the distribution of the local time of pseudoprocesses related to higher-order 
heat-type equations; see Beghin and Orsingher (2005). 

The time- fractional derivative appearing in (1.1) must be understood in 
the sense of Dzerbayshan-Caputo, that is 



■ / 7 r-n 7:. — (x,s)ds, ior m — 1 < v Km, 



— — (xA), for z> = m, 

where m — 1 = [i^J . 

Considering the derivative in the sense of Dzerbayshan-Caputo permits 
us to study initial value problems for (1.1) with initial data represented by 
derivatives of integer order; on this topic, consult Mainardi (1996). 

We assume, in particular, the following initial condition: 

(1.2) u{x,0) = 6{x) forO<i/<l, 
and 

(1.3) |''^f'°ir_'^i'^^' forl<^<2. 



ut{x,0) = 0, 

The general solution to equation (1.1) subject to (1.2) or (1.3) is well 
known [see Podlubny (1999), formula (4.22), page 142] and reads 

i_ - (-Ix|/(AW^))^ 
"^ '^ 2\t^l^ ^^k\T{-uk/2 + l-iy/2) 

(1.4) 

1 



where Wa^p in (1.4) denotes the so-called Wright function, whose general 
form is 

00 J. 

(1.5) ^..(x)=E ^,^J^^^^ , a>-l,/3>0,x.M. 

Some properties of the Wright function are investigated in Mainardi and 
Tomirotti (1998) and in Gorenflo, Mainardi and Srivastava (1998). Initial 
value problems (as well as problems on half-lines with boundary conditions) 
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for equations like (1.1) are extensively treated and solved in Mainardi (1994, 
1995a, 1995b), Gorenflo and Mainardi (1997) and Buckwar and Lucliko 
(1998). 

It has been proved also that Uy is nonnegative and integrates to one for 
all < < 2; see, for example, Orsingher and Beghin (2004). 

We present here some alternative forms of the solution Uy of (1.1), either 
as integral functions like 

or in terms of stable densities 

p„(x;7,7?) = — / e-^^^exp{-r/|/3|°e-^'^^/2/3/|/3||^^^ ^^^^ 
27r ./_^ 



as 



Uu{x,t) 



lp2/u(\x\;lii^-l),X^/''t^, l<i^<2. 

In Orsingher and Beghin (2004), we proved that in the special case v 
the solution (1.4) coincides with the distribution of the process 

(1.6) Zi{t) = Bi{\B2{t)\), t>0, 

called the iterated Brownian motion, which consists of a Brownian motion 
Bi whose "time" is an independent reflecting Brownian motion. 

In Beghin and Orsingher (2003) we have generalized this result to the case 
where u = ^, n S N. In this case, for = 1/2, the solution (1.4) coincides 
with the distribution of the process 

(1.7) Jy^{t) = Bi^J[Gj{t)y n>l,t>0, 

where the vector process (Gi(t), . . . ,Gn-i{t)) has the following joint distri- 
bution: 

(1.8) 

Wj > 0,1 < j < n — 1, 

for n>2. 

In (1.7) the role of "time" is played by the product of independent, 
positive-valued r.v.s, which cannot be identified with well-known distribu- 
tions as in the special case (1.6). 
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In the special case n = 2, we note that Ji/2(t) = Zi(t), because (1.8) be- 
comes the distribution of a reflecting Brownian motion. 

We are now able to prove a much stronger result for the case u = n ^'N, 
and for = 2^/^"~^, which has a number of interesting consequences. We 
will show below that (1.4) for v = can be written down as 

ni/2n (j;, t) = 2" / • • • / — -== -== j=- dzi--- dzn 

Jo Jo v27rzi V2vrz2 \/27rt 

(1.9) 

and this coincides with the distribution of 

(1.10) in{t)=Bi{\B2mi---{\Bn+im---)\)\), t>o, 

where the Bj's are independent Brownian motions. 

The iterated Brownian motion Ii{t) = Bi{\B2{t)\) has been actively in- 
vestigated and many of its properties have been obtained by Khoshnevisan 
and Lewis (1996), Burdzy and San Martin (1995), Allouba (2002). 

The connection between fractional generators of order 1/2 and the iter- 
ated Brownian motion Ii{t) has been studied in Allouba and Zheng (2001) 
and Baeumer, Meerschaert and Nane (2007). This connection was obtained 
in Orsingher and Beghin (2004) as a particular case of the analysis of the 
fractional telegraph equation. 

The identity 



(1.11) 



-xV(22i) p-zf/{2z2) p-zlim 

■ dz\ • • • dZn 







^/2-KZi ^J2^TZ2 y/2TTi 

1 _ - (-2|x|/(2t)V2"+^)fc 



(2i)i/2"+i fc!r(-A;/2"+i + 1 - l/2"+i) 



shows that there is a deep connection between Wright functions and Gaus- 
sian distributions. 

For the n-times iterated Brownian motion X„(t), t > 0, we obtain the 
distributions of the maximum and the sojourn time (together with the ex- 
pression of moments) and we work out in detail an explicit form of them for 
the case of the classical iterated Brownian motion Ii{t), t>0. 

We note that X„(i) converges in distribution, for n — )■ +oo, to a Gauss- 
Laplace (or bilateral exponential) random variable, independent from t > 0. 

In Orsingher and Beghin (2004) we have seen that for the fractional 
telegraph-type equation 

{du d^/'^u 2^^^ 
u{x, 0) = 5{x), 
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the general solution coincides with the distribution of the telegraph process 
T whose time is an independent reflecting Brownian motion 

(1.13) W{t)=T{\B{t)\), t>0. 

We remark that process (1.13) converges to (1.7) in the Kac sense (i.e., 



,2 



for A, c — )• oo, in such a way that ^ — )■ 1). Related interpretations of the 
solutions to 



u 



(1.14) <; Qi2. at- dx^' 

[u{x,0) = 6{x) 

are discussed in Beghin and Orsingher (2003) and Orsingher and Beghin 
(2004). Generalized forms of the fractional telegraph equation (1.14) and of 
its solutions can be found in Saxena, Mathai and Haubold (2006). 

We obtain here various types of relationships between the solutions for 
different values of z/. The first one we present is the following: 

1 

(1.15) u,ix,t) = ^ e-''/^^'^U2Ax,z)dz 



(valid for any < u < 1), where U2u is the solution of (1.1) with order 2u 
instead of u. Formula (1.15) leads, for v = to the n-times iterated Brow- 
nian motion defined in (1.10), since it permits us to obtain, in an alternative 
way, the relationship (1.9). 

In the general case, (1-15) shows that the process related to the equation 
(1.1) of order v can be interpreted as the composition of a process governed 
by the same equation, but with order 2u, with a Gaussian-distributed time. 
We also derive the analogous relationship 

(1.16) u,{x,t) = / e-"'/(^^"'%2.(u;,t)dw, 

Jo yjAirXw 

where 

(1.17) n,.(«.,t) = {^--K^)' ->0; 

Here the roles of space and time are interchanged with respect to (1.15). 
Therefore from (1.16) a further interpretation of the solution emerges, be- 
cause it coincides with the density of the process 

B{%{t)), t>0, 

where i? is a Brownian motion and Tu{t) is a process independent from B 
with a distribution for each t given in (1.17). 

A relationship similar to (1.15) and connecting with Umv is established 
(by applying the multiplication formula of Gamma function) for m > 3 and 
< 1/ < 2/m. 
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Substantially different situations are encountered for the special cases 
1/ = ^ , = I and 1/ = I . In particular for = | we show that the solution to 
(1.1) possesses the following simple form: 

where Ai(x) is the Airy function. The latter emerges as a solution to third- 
order heat-type equations of the form 

du d^u „ ^ 

By using again the relationship (1.15) we get, for the case = ^, the 
following result: 

This suggests that we should interpret Uij^ as the distribution of 

J^/^{t) = A{\B{t)\), t>0, 

where A is a process whose one-dimensional distribution is given in (1.18), 
which coincides with the symmetric stable process of order 1/3. 

Similar relationships seem not to hold for the solutions to fractional equa- 
tions of order u = z:,'n > 3, because the fundamental solutions to 



(1.20) 



du d^u 



dt " 



Cn = ±1, are sign-varying functions on the whole x-axis (while, for n = 3, 
only on the negative half-line), as shown in detail in Lachal (2003). Therefore 
they cannot be used to construct the functions Uy emerging from (1.1), 
which, for < < 2, are nonnegative and integrate to one. We note that 
the solutions to (1.20) themselves have been represented as distributions 
of compositions of artificial processes, which do not display a probabilistic 
structure [see Funaki (1979), Hochberg and Orsingher (1996), Benanchour, 
Roynette and Vallois (1999)]. 

Finally the previous results permit us to establish connections between 
the solutions U2/3n and U2j^n-i. Moreover the explicit form (1.18) of U2/3 
suggests that we should interpret them as distributions of processes similar 
to the n-times iterated Brownian motion, but with the role of B replaced 
by A and the time represented by nested products of the random variables 
Gj defined in (1.7). 
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2. Iterated Brownian motions generated by fractional equations. In this 
section we examine in detail various relationships between solutions to dif- 
fusion equations like (1.1) and processes involving Brownian motion. All 
results of this section refer to equations of order < < 1 . 

We start with the following general theorem: 

Theorem 2.1. The solution to 



u 



(2.1) I 'd¥~^^'d3?' xER,t>0, 
L u{x, 0) = (5(x), 

for Q <v < 1, can he represented as 

(2.2) u,{x, t) = ^ e-''/^^'K2u{x, z) dz 

VTTt Jo 



where U2u is the solution to 



X 



u 



1 



(2.3) <^ Qz2v ^ dx^' forO<u< ^ 

I u{x, 0) = S{x), 



or 



(2.4) 



u{x,Q) = 5{x), for ^<v<i. 
, ut{x,{)) = 0, 



Proof. By applying the duplication formula of the Gamma function we 
have that 

(2.5) r + 1 - - = VvF2^('^+^' ^ y ^ )> 



2 2 y r(l/2(l -z^(A; + l)))' 

By plugging (2.5) into (1.4) we get that 

1 - (-|x|/(At-/^))^r(l/2(l-Kfc + l))) 
' ' 2AW2 ^ ik!0F2-('=+i)r(l - u{k + 1)) 



V^2'^+iAW2f^ fc!2'^^r(l-i/(fc + l)) 

00 

g-^^-./2-l/2 



0F2'^+iAt^/2 Jo 



°° 1 1^1 

^ ^ fc!r(l - z.(A: + 1)) V" A2-(u;t)-/2 J "^"^ 
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[in view of (1.4) with suitable arrangements] 

^l^t^l^ Jo \ J zuy , J 

Jo 



[2^/tw = z] 
1 



oo 

2 / 



e 



'T^t Jo 

and this concludes the proof. 

An alternative proof of the relationship (2.2) is based on the Fourier 
transforms, since for Ui, the following result is known: 

e''^''u^{x,t) dx = E^^ii-P^Xh"), 



oo 



where Ey^i{z) = X^fcLo r(fct+i) Mittag-Leffler function. Taking the 

Fourier transform of (2.2) we get that 



-hoo 

e 

oo 



'^^l-^J'^ e-'"'/^^'^U2uix,w) dw^ dx 



1 



oo 

2 



e-^'/(^'^E2,i{-P'XWndw 

'irt Jo 

: > -V^^ ^ / ■^^w^'^" dw 

^^T{2ky + l)Jo V^t 

: [for w = 2Vtz] 

r('9t;y 4-1 O,/:^ \ 9 



A;=0 

oo 



i-p^x^)^ 


(2Vt)2fc'^+i 


T{2ku + 1] 


1 2V7rt 


{-p^x^f 


(2^*)^^''+^ 


T{2kiy + l] 


1 2^/7rt 


(-/32aV) 


fc r+oo 

= / e^' 


V{kv + l) 


J —oo 



,.0 "V... -X^fc) 

oo 

E 

fc=0 

Remark 2.1. In the special case where ^, formula (2.2) yields 

Uy2{x,t) = ^ e-^/(^*) " dz 

Virt Jo VAttX^z 



□ 
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2, 



(2.6) = [2\'z = y] 

Particularly interesting is the case where 2(2A^)^ = 1, that is, when = 
2-3/2^ because (2.6) reduces to 

(2.7) u,M^,t> = 2l ^ ^ ,y, 

which permits us to conclude that, in this case, the solution coincides with 
the probability density of the iterated Brownian motion (1.6). 



Remark 2.2. If we generalize our analysis to the n-dimensional case and 
1 

2' 



take 1^ = h, we can show that the process related to a fractional equation of 



the form 

with initial condition 

n 

o)=n'^(^*^)' 

k=l 

has components represented by iterated Brownian motions with a common 
random time. In other words, the solution to (2.8) coincides with the distri- 
bution of the vector process 

Bimm, 

t>0, 

Bn{\Bm, 

where B^, k = l,...,n, are mutually independent Brownian motions and 
also independent from B. 

To check this result we evaluate the Fourier transform of the solution to 
(2.8) as follows: 

f + OO /• + 00 

e'^^"i+-+'^"""ni/2(a;i, . . . ,Xn,t) dxi ■ ■ ■ dxn 



oo 



(2.9) =^iAi(^-A^i^/'(^EA 
2 



10 E. ORSINGHER AND L. BEGHIN 

From (2.9) we get the inverse Fourier transform in the fohowing form: 

= 2 1 — — I I — , dw. 

Jo ^/M2W) V2^ 



The main difference with respect to the case of the usual multivariate heat 
equation is that the components of the iterated Brownian motions are no 
longer independent because they are related to each other by the common 
random time B (with infinitesimal variance 

We pass now to our second theorem, which is related to the case v = i^, 
n G N. 

Theorem 2.2. For = ^, A = 2^/^'"^^'"^ the solution to equation (1.1) 
under the initial condition (1.2) can be written as 

ui/2n(a;,t) = 2" / ••• / — -== -== -^=^dzi---dzn. 

Jo Jo V^'^zi \/2txz2 V27rt 

(2.10) 

Proof. In view of the duplication formula for the Gamma function we 
can write 



(2.11) 



2«+i 2""'"-'- 



/3oi/2"+fc/2" r(l - A;/2" - 1/2") 
^ r(l/2- A:/2"+i -l/2"+i) 



so that the first member of (1.4) becomes, for 1/ = ^ and A = 21/2'"+'' -i^ 
Mi/2n(x,t) 

J 2\x\ y 

)l/2"+iZ^\^ (2t)V2"+ij 



(2i)i/2"+i Z^\^ (2t)V2"+iy fc!r(l - A;/2"+i - l/2"+i) 

1 ^/ 2\x\ Yf.-e-^w^'/'-^'-'/'-'^'-'/'dw 



1 y./ 2\X\ \ 

l/2"+iZ^\^ (2t)l/2"+iJ 



(2t)i/2"+i (2t)i/2"+iy A;!V^2(fc+i)/2"r(l-A;/2"-l/2") 

1 

(2t)i/2"+i 
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coo -l/2"+i-A:/2"+i-l/2 , roo _^„„ -l/2"-fc/2"-l/2 , 

Jq e ^tM;^ dwi Jq e dw2 



1 



(2t)i/2"+i 

2|x| 



oo 

X 

A;=0 



E 



(2i)V2"+i 



fc!(V^)"2('^+^)S^=o ^/2""'r(l/2 - k/2) 
At this point we can use the reflection formula for the Gamma function 



(2.13) 



2 2 J sin{7r/2(l-fc)}r((l + /c)/2) 

TT 



cosA;7r/2r((l + fc)/2) 



and this shows that only even terms of (2.12) must be retained. We can 
therefore write that 

1 



(2t)i/2"+^ 

2\x\ 



oo 

X 

k=0 

n 



E 



(2t)i/2"+i 

xcos^r(i±^)[fe!(V^)"7r2('=+i)S;^"oV2"-pi 



E 



(2t)i/2"+iz_.^ (2t)i/2 



(2.14) X 



2|x| 

1+1 I 

(r • • • r ^- n -^r ^-(^+i)/2"^---i/2 



/fc!(V^)"+i2Si-oi/2"-2'=Si-oi/2"-r(fe/2) 
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(2t)i/2 
oo 



r=0 



2\x\ 



2r 



(2t)i/2"^ 



r! 



( V^)-+i2^i=o ^2"-^- 22-E;.o 1/2" ^ 
2 



(2t)l/2"+i22(l-l/2"+i)(^)n+l 

{-ly ( x2 



E 

r=0 



^2-1/2' 



! V 



(2- 



~2(2-l/2")jr 



oo fOO 



e-E"..-.-Q^-(2^+l)/2"+--l/2^^^. 



By considering that 



(-ir 



r! 



a; 



r,\ 1/2" n 

2 \ ' T-r -l/2"+i 



r=0 



we can write (2.14) as follows: 

1 



-1/2" 



^l/2"(x,t) = ^^^ 



' 21-1/2" (^) 



n+1 



00 fOO 



,2/22(2/4)1/2" n".i«': 



i=i 



-1/2' 



n+2-j_ 



-1/2 



-1/2 



n + l—i 



In order to calculate the integrals let us write 



2(2-it)V2"n 



w 



1/2 



n-l-l— j 



Zl 



i=i 



so that 



and 



" \ 2 



zi2-i(2-H)-V2 

T-rn-l l/2"+l-i 

ni=i 

„_l/„_l,x_l/2" X 2 



PROCESSES WITH RANDOMLY VARYING TIME 

Therefore we get 

1 

Ui/2r^{x,t) - ^2^^1/2"+! 21-1/2" (^)n+l 



OO fOO 



^0 



n-1 



n 



X 2zi 



2-i(2-it)-i/2"\2 



1-rn-l l/2"+i-J 

n,=l 



J=l J 

^i2-i(2-H)-i/2"\ -1-1/2 



ni=i 

n + 1 



(2.15) 



dzi dwi ■ ■ ■ dwn 



(2t)i/2"+^ 21-1/2" (^)«+i 



zi Jo 

X e- S;=i' JJ ^-1/2"+!-^ -1/2 _ _ 

i=i 

Now we make the similar substitution 

n-l 

i=i 

so that we get again 

■222-i(2-it)-i/2"-^2 



W^n-1 

and 



nn-2 1/2"-^ 
.=1 "^i 



dw„_i = 2z2 dz2 



2-i(2-it)-i/2"-^y 

T-r"-2 l/2"-J 

n,=i 



^J = l J 

In view of these substitutions, formula (2.15) is transformed into 

_ V2(2-it)-i/2"+^^(2-it)-i/2" 
(2t)i/2"+^ 21-1/2" (0F)«+i 
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roo ^-x^/{2zi) roo ^-zf/(2z2) 
(2.16) X / — dzi / — dz2 

Jo \fz\ Jo 



By similar transformations, after (n — 3) additional steps, we arrive at 

_ V2^(2-lt)-l/2"+^-l/2" — -1/23 
^1/2" = (2t)V2"+i2l-l/2"(^)n+l 



dzi / — — dz 



:2 • • 



/^^l Jo V^2 

'^-4_./(2^[(2-*)V2.y^])^_„,^-l/22-l/2^^^^ 

By means of the position 

2{2-hf l^w\'^ = zn 

we get that 

w, = {zn2-\2-hr'l'f 

and 

dwi = 2zndzn{2~^{2~^t)~^''^f. 
We arrive at the final expression 

/2^(-2-li)-l/2"+^-l/2" 1/23-1/22 



(2t)l/2"+i 21-1/2" (^)n+l 

^oo g-x2/(22i) r-oo ^-zl/{2z2) 

- dz^ / — dzo 



V^i Jo V^2 



V^' 

2" 



2"/2+l/2(^)n+l^ 

00g-x2/(2.i) .OO g-22_,/(2.„) 

-dz^-- - / — e ''"'^^^'dzr 



/o V^i 
which coincides with (2.10). □ 
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Remark 2.3. It is well known that the Laplace-Fourier transform of the 
solution to (1.1) with initial conditions (1.2) or (1.3) is equal, for < < 2, 
to 



(2.17) 



' + 00 



+ 00 



e*^^nj,(x,t) dx 



s>0,/3g 



We check that the Laplace-Fourier transform of (2.10) reduces to (2.17) 
for 1/ = ^ and = 2^2-2: 



= 2" 



= 2' 



ni/2" {x, t) dx 



e'^"" dx 



-oo 
oo 



oo g-a;2/(2zi) 



\/27rZi 



dzi ■ ■ ■ 



oo g-z2/(2t) 



-^?/(2^2) 



\/27ri^ 



r=0 

oo 



r- ]^ |.oo ^^-zl/{2z2) 

Y I r! 



3 4- 



dz\ ■ ■ ■ 



V2^ 
dz2 ■ ■ 



dZn 



/27rt 



■ dZr, 



/27rt 



r=0 



2 y r! 

72 6 



■oo , ^-z^/(2z3) 



2 

(iz2 • • 



r=0 



'2 ^2^ 

p2\r-^ 2^/2-12^/4-1 



■oo p-4/(2t) 



/27ft 



■ dZr 



■°° . e~^3/(224) 



- + - r 

2 2 



r 1 
4 + 2 



r-/4' 



/27rt 



■ dZr 



r=0 



/32\ 1 2^'/2-l2^'/4-l . . . 2^/2""^-l 



2 / r! 



(\/^)-i 
1~ 



r 

22 + 2 



r 1 

2"-i + 2 



-4/(2t) 



r=0 



V27rt 

^2\ r- 2'"/2+»'/4+-+r/2"- 



.r/2" 



2 / r! 
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By applying the duplication formula we get that 
(2.18) = ^2.-'^V?2'-'/^iM . . . ^2-'/-- 



r(r/2) r(r/22) ^ r(r/2'^ 



r/2"-i ^(r) 



r(r/2" 
and thus 

-oo 



oo 



r=0 

oo 



f^'^\ J_2r-/2+r/4H \-r/2^-n2n-r-r/2 r/2"-i^r/2" ^X0_ 

r! r(r/2") 



QnV^/ \ _!_2r-/2+r/4H \-r /2" -rirfU-r-r /2 r/2"-i^r/2 

2 y r! 

^2\ r 2r-/2"-r^r/2" 



r=0 ^ ^ 

oo 



r=0 

E- 



r/2"r(r/2") 



22-1/2" y r(r/2" + l) 



1/2'M 1-22-1/2" 



By taking the Laplace transform of (2.19) we get 

r'+oo / ^2^1/2" \ ^l/2"-l22-l/2" 



e ''*^l/2«,l(^-22_l/2" j ^2 + 22-1/2" sl/2"' 

which coincides with (2.17), for = ^ and = 2^/2" -2^ 

The form (2.10) of the solution iii/2" shows that it coincides with the 
distribution of the n-times iterated Brownian motion defined in (1.10). 

Another representation of the solution to the fractional equation (1.1) can 
be inferred from the following result: 

Theorem 2.3. The solution Ui^{x,t) = Uu to the initial value problem 
(2.1), for < I' < 1, can he written as 

(2.20) u^{x,t) = / e-^'/(^"'^%2^(w,t)d«;, 

Jo v47ri(;A 
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where 

(2.21) ^2.Kt) = {'"^^("'*)' {"'""l'' 

[0, for w <\J 

and U2u is the solution of (2.3) or (2.4)- 

Proof. We first note that for tlie solutions to (2.3) or (2.4) the following 
result holds: 



/•oo 

(2.22) L{x,s)= / e-'^U2u{x,t)dt 

Jo 



2A 



as can be obtained by taking the Laplace transform of = A^j|. The 
solution to the corresponding equation 

coincides with the solution to 



dL 
dx 



dL 
dx 



s 



A2 ' 



L(s,0+) = L(s,0-), 

and easily yields (2.22); see also (3.3) of Orsingher and Beghin (2004). There- 
fore, by taking the Laplace transform of (2.20), we get 

^^g-x2/(4^A) L r e-''u2u{w,t)dt\dw 

v4vri«A [Jo J 



V4:TrwX 2A 
[2w = z] 



2A Jo y/2TTzX 



1_ g-|x|WVA 



e 



x^/(2zX)^-s-^/Xz/2 



2A 

and this coincides with the Laplace transform of u^{x,t). □ 

Remark 2.4. Formula (2.20) suggests that we should represent the so- 
lution of (2.1) as the distribution of the process 

(2.23) B{T2u{t)), t>0, 
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where B is a. Brownian motion with infinitesimal variance 2 A and 72u{t), 
t > 0, is a process, independent from B, with law equal to (2.21). 

It is straightforward that, for u = 1/2, the process (2.23) coincides with 
the iterated Brownian motion Ii; see (1.6). 

By comparing the relationship (2.20) with (2.2) we note also that, in the 
composition of processes, Brownian motion plays in the second case the role 
of "time," while in the first one it represents "space." 

3. On moments and functionals of the iterated Brownian motion. Some 
properties of the classical iterated Brownian motion have been obtained 
by several authors and include the law of iterated logarithm [Burdzy and 
San Martin (1995)] and the modulus of continuity [Khoshnevisan and Lewis 
(1996)]. Applications of the iterated Brownian motion to diffusion in cracks 
are dealt with in De Blassie (2004). 

We start by presenting the distribution of the maximum of the n-times 
iterated Brownian motion and, in an explicit form, for the usual iterated 
Brownian motion. 

Theorem 3.1. For the n-times iterated Brownian motion 

Xn{t) = Bi{\B2{\B^{- ■ ■ {\Bn+im ■ ■ ODD, t > 0, 

where Bj, j = 1, . . . ,n + 1, are independent Brownian motions, we have for 
/3 > that 



Pr< max X„(s) G dB 

' 0<s<t 



(3.1) 



I-+00 /"+00 r ■\ 

2 ■■■ Pr{5i(yi)Gd/3}Pr<^ max |S2(zi)| G dyi ^ 



X Pr<^ max \B-i{z2)\ ^ dy2 \ ■ ■ -Vil max |B„+i(2;„)| G (iy^ 

0<Z2<?/3 I I 0<^„<t 



Proof. For Xi{t) = Bi{\B2{t)\) we can write that 



Pr< max Xi{s) G dB 
' 0<s<t 

= Pr| max Bi{z) £ dfi) 

[0<z<maxo<u,<t |-B2(w)| J 



(3.2) =E{Vt{ max Bi(z) e dP max 1^2 (w) 

[ [o<2<maxo<„<t|B2(«')| 0<«><t 

= / Pr| max Bi(z) G dlsX Pr| max \B2(w)\ G dy\ 
Jo [o<z<y J [o<w<t' J 
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= 2 [ Pr{Bi(y) G d/3}Pr| max \B2{w)\ e dy\ . 

Jo [0<w<t J 

For = i3i(|X„,„i(t)|), n > 1, we have analogously that 
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(3.3) 



Pr< max € dB 

' 0<s<t 



= 2 Pr{5i(y) G d/3}Pr| max |X„_iH| G dy| 

and, by induction, we obtain (3.1). □ 

Remark 3.1. In the case n = 1 we can give an explicit expression for 
(3.1) as follows: 



Pr< max Xifs) G dB 
' Q<s<t 



■ 2dp 



+00 -/32/(2t«) 



\/2ttw 

( +00 



(3,4) 



, fc=— oo 



[l + 2k)- 



-«>V(2<)(l+2fc)2 

y/2Trt 

-u>2/(2t)(l-2A:)2 



'2-Kt 



+ {l-2k) 

:2g(-l)'=[pr{x,(^^).d/5KPr^^i 



fc=— oo 



dw 



t 



;i-2ifc)2'^^^ 



2ky 



where Ui/2{x,t) is given in (2.7) and in the first step we applied the well- 
known result for the maximal distribution of the absolute value of Brownian 
motion [see Shorack and Wellner (1986), page 34]. The last term of (3.4) 
shows that the distribution of the maximum of the iterated Brownian motion 
can be expressed in terms of its probability law Ui/2 — ^1/2(^1^)) the 
case of the classical Brownian motion. 

In principle we could write explicitly the distribution of the maximum of 
Init) in terms of but this produces a sum of 2" terms, each of which 

has a very entangled structure. 
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On the basis of the same principles it is possible to write down the 
distribution of the sojourn time on the positive half-line of the process 
= t > 0, n > 1, defined as 

(3-5) ^t= '^{z:Biiz)>0}dz. 

Jo 

This random variable takes values in [0, +oo), because during the interval 
[0,t) the process (which plays the role of time for Bi) can span the 

whole positive real axes. 

Theorem 3.2. For the process X„(t), t > 0, the distribution ofVt reads 

PrjE/ G dsj = ds / Pr< max \Zn-i{w)\ & dz > , 

^ Js ^Js{z-s) \o<u,<t^ ^' j' 



(3.6) 



< s < oo. 



Proof. The definition of Tt given in (3.5) implies that 
Pr{ri G ds] 

fmaxo<u,<t \Tn-i{w)\ 



(3.7) =E{Vt< 



l{z:Bi{^)>0} dz 



G ds 



max \Tn^i{w) 
0<w<t 



= Pr{r^ G ds}Pr|^maxjX„_i(w)| G(iz|. 

By inserting the arc-sine law in (3.7) we get (3.6). 
We can check that (3.6) integrates to one 

Pr{rt G ds} 

r+oo r+oo ^ f 

= / ds Pr< max \Zn~i(w)\ ^ dz > 

Jo Js TT^S{Z-S) {0<w<t J 

= / Pr< max \In-i{w)\ £ dz > / . = 1. ^ 

Jo lo<«-<t' " ^' j Jo 7r^s{z-s) □ 



Remark 3.2. For the iterated Brownian motion Ii{t) = Bi{\B2{t)\) the 
distribution of Tt can be written explicitly as follows: 

M^t e ds} 

"+°° dz 



ds 



TTy/s{z - S) 
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(3.8) X J2 (-1) + + 7t=. (1 - 

k= — oo 

J +2° C /.+00 -22/(2t)(l+2fc)2 



+00 g-z2/(2t)(l-2fc)2 



+ (1-2A:)/ ^ dz^. 

V2; - s J 

By the transformation z = s{l + x^) the integrals in (3.8) are converted 

z2k 
2t 

+00 -s^A(l+x^)^ 



[for ^ = 11^] into 



2s / ^ dx 

\/s 



Jo 



+ 00 



where, in the last step, we have applied formula 3.469.1 of Gradshteyn and 
Rhyzik (1994) and K^/^ix) = -^[I_y^{x) - h/^ix)] [by formula 8.485 of 

Gradshteyn and Rhyzik (1994)]. By we denote the Bessel function of 

imaginary argument of order u, that is, Iu{x) = X]/t^ fc!r(fc+i^+i) ' Therefore 
we get 

Pr{rt G ds} 

/ ,-,kfn , ni.^.~sVm(l+2k)^r. fs'^{l + 2k) 



(-l)'=|(l + 2fc)e-^'/(^*)(^+2A:)2i^^ 



k=—oo 

+ (1 - 2fe)e-V(^*)(i-2^)^i^,/4(^^^i^) }. 

We now derive the explicit form of the moments of even order of X„(t). 

Theorem 3.3. For the process Xn{t), t > 0, the moments of order 2k 
are given by 

(2kV2'^ f°° , foo -x^/{2zi) pco -zl_J{2t) 

k\ 2*- Jq Jq y/2TTZi Jo V2vrt 



(3.9) 



2^/^" (2fc)! 
22fc r(A;/2" + l 
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Proof. The first expression in (3.9) can be proved by observing that, 
for n > 1, 

El'Jit) = E[Bl''{\B2m{---\Br.^i{t)\- 





{2k)\ 


1 


(3.10) 


k\ 








{2k)\ 


1 




k\ 






{2k)\ 


1 




k\ 


2k 



E\B,{\B,{---\Br.+,{t)\-.-)\t 

r+oo 

2 / x'=Pr{52(|S3(---|5n+i(i)|---)l)edx} 

JO 

r+oo 

2 / x''Vv{Za^i{t)edx], 
Jo 



which coincides with the second hne of (3.9). By performing the integrations 
in (3.10) we get the exphcit expression of the moments of order 2k: 

r{k/2 + l/2)r(A:/22 + 1/2) • • • r(A;/2" + i/2)2'^/2+-+fc/2"tfc/2" 



X 



2"V7r'' 

(2fc)!2" 



k\2>' 

(3.11) = [by (2.18)] 



r(A:/2") 2^^^ A:! 



. ^fc/2"2n-2fc+fc/2" 



{2ky. 



kT{k/2'^)' □ 

Remark 3.3. For n = formula (3.9) coincides with the moments EB'^^{t), 
which is as it should be, since Xo(t) = B{t). 

For n= 1, the moments of the iterated Brownian motion Zi{t) = Bi{\B2{t)\) 
can be evaluated directly as follows: 

Elf\t) = EBlWB2m 
{2k)\ 1 



k\ 2^ 



E\B2{t)\' 



kl 2^ V2^ 

^k/2 



2^1^ {2k)\ 



22fc r(A;/2 + l) 

which coincides with (3.9) for n = 1. 

For any n > 1 and A; = 1 , we obtain the explicit form of the variance 

2l/2"^l/2" 

varX„(t) = — — — ; -, 

^' 2r(l/2" + l)' 
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while, for n = 0, it is varXo(t) = t, as expected. 

Remark 3.4. For all t > 0, the sequence Tn{t) converges in distribution, 
for n— 7- oo, to the Gauss-Laplace exponential random variable and its den- 
sity is independent from t. From (1.11) we get that 

(3.12) lim ui/2«(x,t) = e"2|'^l, i>0,xGM. 

n— >oo ' 

By working on the Fourier transform (2.19) of ^^^2" we have the following 
alternative proof: 

Formula (3.13) coincides with the characteristic function of (3.12). Loosely 
speaking, this shows that the composition of infinite Brownian motions pro- 
duces the bilateral exponential distribution. 

In view of (1.9) we have also the identity 



(3.14) 



l-oo r-oo g-x2/(22i) ^-zl/(2z.2) g-4/(2i) 

lim 2" / • • • / — dzi--- dzn 

Jo Jo v2vrzi V27rz2 V27ri 



which is a rather striking result. Furthermore, if we assume that 
the fractional equation (1.1) is converted into 



1 



22 dx^ ' 

subject to 

u{x, 0) = 5{x), 
which is satisfied by (3.12) for all x / 0. 

Remark 3.5. For the random process 

(3.15) r(|i?2(|i?3(---|^n+l(t)|---)l)l), t>0, 

where T is a telegraph process (with parameters A and c) independent from 
the Brownian motions = 2, . . . , n -|- 1, we have a similar result. The 

distribution Ui/2"- of (3.15) is a solution to 

92/2" 5^/2" 2^^"" 

u{x, 0) = 5{x), 
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and its characteristic function is equal to 

"+00 



(3.16) 



e*^"" til/an 


1 




~ 2 







A 



- c2/32 

A 



VA2 - c2/3' 



^l/2M 



where r^i = — A + y^A^ — c^/S^ and r/2 = — A — y^A^ — c^/S^ [see Orsingher and 
Beghin (2004), formula (2.7), for a = 1/2"]. 
For n — )• oo we get from (3.16) that 



(3.17) 



lim 

n— >oo 



Ui/2ri{x,t) dx = — 



1 + 2A 



+ 2A + c2/32^ 



which is the characteristic function of the bilateral exponential random vari- 
able, with density 



(3.18) 



/(^) = Vl_t2A^_|,|^I+2A/c^ 



X G . 



Clearly, for A = and c = 1/2, (3.18) reduces to (3.12) and (3.17) coincides 
with (3.13). 



4. The explicit solution of the fractional diffusion equation for = 1/3, 
f = 2/3 and = 4/3. In some special cases it is possible to present the 
solutions of the fractional equations (1.1) in a more attractive fashion. This 
is the case for v = ^. The explicit form of U2/3{x,t) is given in the next 
theorem, in terms of Airy functions. 

By combining this result with the relationship given in Theorem 2.1, 
Ui/-^{x,t) can be represented consequently in an interesting form. 



Theorem 4.1. The solution to 



(4.1) "(9^273" xGM,t>0 

u{x, 0) = 5(x), 



can he represented as 
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where 



Ai{w) 



(4.3) 



vr 



+00 



cos I aw + ^ ] da 



w 



1/2 



is the Airy function and ly denotes the Bessel function of imaginary argu- 
ment of order v. 

Proof. From (1.4) we readily have that 

{-\x\/{\t^/^)f 



U2/3{x,t) 



1 (-F 
vtV3 k\T(l 



(4.4) 



2AtV3 A;!r(l - (A; + l)/3) 

fc=0 



1 



27rAtV3 



E 

fc=0 



(-|x|/(Ai^/3))fcp((^ ^ l)/3) sin(^(fc + i)/3) 



A;! 



By direct inspection the following identity is proven to hold: 
(4.5) 



sm 



3 ' ' 3 
and, by inserting this into (4.4), we get that 

(4.6) ,3(.. t) = ' E + '"'P-'^ + . 

V ; 2/3\ ) 2^AtV3Z^ A;! 

fc=0 

We note that, from (4.3), for all \w\ < 00, 



Ai{w) 



w 



1/2 



3 

^1/2 



/2wrV2 



/2u;3/2 



■^/^/2x2fc-l/3 



.fc=0 



3 
1 



fc!r(/c- 1/3 + 1) 



(4.7) 



E 

A;=0 
3k 



W 



3/2\ 2fc+l/3 



1 



fc!r(A: + 1/3 + 1) 

,,3fc 



E- 



w 



32k+2/3klT{k + 2/3) ^ 32^-+4/3 A;!r(A: + 4/3) 

k — A? — 



E( ^ ) 



sin(27r(fc + l)/3) 



37/6 zL. ^32/3 J r{ik + 2)/3)r((fc + 3)/3) ' 
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The last step can be justified by taking k = 3m, 3m + 1 and 3m + 2. While 
for k = 3m + 2 the last term in (4.7) is equal to zero, in the other two cases 
the two series are obtained. 

The triplication formula of the Gamma function [see Lebedev (1972), page 
14], that is, 

(4.8) riz)r(z + l)r(z + l) = :J^ri3z), 

k+l 



3 / 33^-1/2 



for z = ^ yields 

(4.9) ^(k_ + 2\fk + 3\_ 2vr r(A; + l) 



3 J \ 3 J 3^'+i/2r((A: + l)/3)' 
From (4.9) we have that 

(4.10) Ai(uO = 5:^V(3V3„)'2H(M^l±iMr('^V 

TT ^-^ kl V 3 / 

and (4.2) easily follows by comparing (4.10) and (4.6). □ 

Remark 4.1. The expression of 1*2/3 {^^ *) obtained in the previous theo- 
rem can be recognized (up to the factor 3/2) as the solution of the third-order 
heat-type equation 



dv . 3 9^ 



(4.11) ^ dt~ ^ 5y3' yeM,t>0, 

v{y,0) = 5{y), 

evaluated at y = Since Ai{y), for y > 0, is positive- valued [see Figure 
1(a)] and the function (4.2) integrates to one (as we show below), U2/s{x,t) 
is a true probability distribution: 

»+oo 

U2/3{x,t)dx 



V 



^ Ai(—^]dx+l -^=Ai{ T^]dx 



Xx/3t \Xx/3tJ J-^oo\</3t V X^3t 



3 /■+°° 

Ai{y)dy = l, 



where the last step follows by noting that J^°° Ai{y)dy = 1/3; see Nikitin 
and Orsingher (2000). 

Therefore we can think of U2/s{x,t) as the probability law of a process 
A{t),t > 0, whose distribution at time t is obtained from the solution v{x,t) 
of equation (4.11), as follows: 

^^2/3(2;,*) = ^v{\x\,t). 
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Fig. 1. The Airy function and the function Ai{\x\) . 

Remark 4.2. For the case = \ the solution uii^{x,t) to (1.1) can be 
written, thanks to the relationship (2.2), as 

1 

(4.12) uy,{x,t) = ^l e-^'/^^'\y3{x,z)dz 



T^t Jo 



1 



32/3 



dz. 



We can represent (4.12) as the distribution of the process 

Jy3{t)=A{\B{t)\), t>0, 

with A and B independent. The results (4.2) and (4.12) show that the 
solutions U2/s{x,t) and Ui/^{x, t) are both unimodal with maximum at x = 0; 
see Figure 1(b). This is in accordance with the general result that, for < 
v <1, the solutions to the fractional equation (1.1) have a unique maximal 
point at X = 0. 

We consider now the case v = 4/3, which is qualitatively different from 
those dealt with so far, because the solutions of fractional equations of order 
1 <u <2 display a substantially different behavior. 

Theorem 4.2. The solution to 



(4.13) 

is given by 

Uii^{x,t) 

(4.14) 



54/3^ 



5*4/3 ' 

u(x, 0) = 6{x), 
I ut(x,0) =0, 



X G M,t > 0, 



1 / 3 \ ^^"^ 



e-'^w-^'^Ai 



\x\ 



A \t 



2/3 



dw. 



28 E. ORSINGHER AND L. BEGHIN 

Proof. From (1.4) we have that 



1 °° / \x\ \ 



2At2/3Z^\ Xt2/3l A;!r(l - 2/3(A; + 1)) 
fc=o 
(4.15) 

1 ^/ \x\ y r(2/3(fc + l))sin(27r(fc + l)/3) 

w2/3 2^1 \+2/3 I h\ 



k=0 ^ ^ 

By means of the duplication formula for the Gamma function we have 
that 

r I i(A: + 1) + 1^ - ^^F2l-2/3(fc+l)^(2/3(fc + 1)) 



r ' 2) r((fc + i)/3) 

and therefore 7x4/3(2;, t) can be rewritten as 
1 



2A^0F2i/3t2/3 



00 



E 

fc=0 



hOO 



X > / e -^i/3{fc+i)+i/2-i 







, ri ^V-^(?)''Y'22(M|±iM,„ 



1 



2A7rV^2i/3t2/3 



°° /■+00 / 1^1 /o\ 2/3 



fc=0 



X > / e-"'«;i/2-2/3 - ^1/3 



A 

sin(27r(A: + l)/3)^/A;+l\ ^ 
" ^l^J'^" 

[by (4.10)] 
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□ 



We can show that 1x4/3(2;, t) dx = 1. Indeed, from (4.14) we have that 



A V t V 3 J 

2 M'^' 



A V t V 3 



Ai[ --( -^/- I 



dxdw 



by the substitution y = ~ — \-\/-^ 



2 /3y/\/2 /T 

v^U; Uv3 

1 2 /3\^^^2~^/^ /■■^°° 



A 

'2/3 p+00 



2/3 



-1/2 



'?3^4 



3-V3 7o 



Remark 4.3. In view of Theorem 2.1 we have the following representa- 
tion for n2/3(x,t), which is alternative to (4.2): 

1 

(4.17) u2/3{x,t) = -y= [ e-^"/^^'K^/,ix,z)dz 



32/3 



2AV^2V3 



+00 



dz 



x\ [ 2 [w\ 

A ViV 3"; 



2/3 



dw. 



By inserting (4.2) into the left-hand side of (4.17) we obtain that 
3 



-.Ai 



\x\ 



2Av^'"VAv^ 

32/3 



r.+oo g-22/(4i) 



dz 



2AV^2i/3 ^^2/3 
[by the substitution s = \/2^twz^] 



(4.18) 



32/3 



2AV^2 



1/3 



^+00 g-^2/(4t) 3^5/3 
^^2/3 25/3t5/6 
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X</3t 



35/3 f+oc / II X r+00 



35/3 ^3/2 / 

— ds 



32/3 



Fr{\B{To)\£ds}Ai 



where 



Fv{\B{To)\£ds} = —-^-^ds, s>0, 



3 5^/2 

27r 1 + s3 

is the McKean law representing the distribution of the position of a Brownian 
motion B at the instant 

r-t 



ro = inf|t>0:l + ^ B(s)(is = o|; 



see McKean (1963). 

By setting y = ^fj— in (4.18) and performing some simphfications we get 

(4.19) Ai{\y\)= Fr{\B{To)\eds}Aii-\y\s), y G M. 



Formula (4.19) shows an interesting property of Airy functions: The value 
of the exponentially decreasing part of 74i(|y|) can be obtained by averaging 
its oscillating component Ai{—\y\s) with the well-known density of |i?(ro)| 
(see Figure 1). 



Remark 4.4. The solution U4^/'^{x,t) can also be expressed in terms of 
3 

2- 



a stable density of order |. Indeed, by using the representation of the stable 



density below 

1 r+00 

(4.20) p,(x;7,r?) = — / e-^^^exp{-r?|/3re-^-^/2/5/l^l}d/3, a/1, 

we know that for aG(l,2), 7^ = 1 and for x > the following series repre- 
sentation holds true: 

(4.21) p.(x;7, 1) = I g;(_,).^^i ^i^(^-(7 + «)/(2«)) r^ + h.\ . 

^ fc=i ■ ^ 

see formula (6.9), page 583 of Feller (1971) (up to some corrections) and 
Lukacs (1969). 
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For a = 4 




2 formula (4.21) reads 




P3/2 



sin{(r + l)2/37r} 
(r + 1)! 




r(|(. + i)) 



If we compare (4.15) with (4.22) we get that 




A different proof of the relationship between stable laws and the solu- 
tions of fractional diffusion equations, based on the inversion of the Fourier 
transform, can be found in Fujita (1990). 

Formula (4.23) proves the nonnegativity of the expression (4.16), as a 
function of x. 

5. Some generalizations of the previous results. In this section we present 
some generalizations of the results of Sections 2 and 4. 

We start by giving a relationship between the solutions Ui, and Umu, Tn>3, 
and obtain some explicit expressions for m = 3. In this case the interpreta- 
tion of U2/3n as the distribution of compositions of different types of pro- 
cesses is possible. Also in this case we encounter processes with a random 
time which possesses a branching structure (depending on n). 

We now state a general result which is alternative to (2.2) and permits 
us to exploit the explicit expression of Uy{x,t). 

Theorem 5.1. The solution to the initial value problem (1.1)-(1.2), for 
< ;/ < 2/3, can he represented as 



(5.1) 




where u^u{x,z) is the solution to 



(5.2) <^ Q^Zu ^ 0^2 

I u{x, 0) = 5{x), 




and 



J dz^'^ dx'^ 
] uix,0) = 6ix), 



2 



x€R,z>0, - 




. ut{x,0) = 0, 
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Proof. In view of the triplication formula (4.8), for z = ^ — ^^^^y^, we 
have that 



-|x|/(AW^))^ 
2AW2 Z^fc!r(l-i/(fc + l)/2) 

k=0 

1 



1 °° 



E 



X 
k=0 



3i-3/2Kfc+i)-i/2r(2/3 - u{k + 1) /2)r(l/3 - u{k + 1) /2) 



A;!r(l-3z^(A; + l)/2) 

V3 



2233/2^'A^W2 
"+00 /■+00 



Jo 



X / / e — ^y;-/2-l/3^-./2-2/3 



X 

fc=0 

^/3 '■+°° '■+°° 



E 



JO 



A(v/3W)3-/2y fc!r(l-3i/(A; + l)/2) 



e-«'"Zy;-i^/2-l/3^-i./2-2/3 



2^33/2^i^/2 

(5.4) X {\/3^wztf''^^U3uix, \^33wzt) dwdz 



27r33/2^t'^/2 



Jo 



g -«;-2y^-i./2-l/3_^-i./2-2/3 



V3 /■+°° /■+°° 
2^ 



JO 



e 



X (w2;)^/^U3y(x, V^wzt)dwdz 
-«,-^^-i/3^-2/3^g^(^^ 3 vW) du- dz, 



which reduces to (5.1), after the change of variables 
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It can be easily checked that, also in this form, the solution integrates to 
one. By using the last expression in (5.4) we get 

f + OO 

Uu{x, t) dx 

oo 

/o /"+00 f+oo r+oo 

= — / e-'"^"u;-i/3^^2/3 / u-i^{x,2,^f^t)dxdwdz 

27r Jo Jo J_oo 

2vr io Jo 

since, by the triplication formula for z = 1/3, it is r(|)r(^) = 2tt/\/^. 

Remark 5.1. By using the previous result it is possible to obtain alter- 
native forms for the solution to the initial value problem for = 1/3 and for 
L' = 2/9. Indeed, in the first case it is 

uy3{x,t) = ^^J J se-^''+'''^/^^'^'Ki{x,sv)dsdv 

(5.5) 

hoo f+OO ^ ^ ^ p-a:2/(4A2{si,)) 



27r\/t Jo Jo 2\y/¥sv 
The relationship (5.5) shows that ui/3 can be interpreted as the distri- 
bution of a Brownian motion (with infinitesimal variance 2A^) at a random 
time Gi{t) ■ G2{t), that is, 

(5.6) Jy,it) = B[Giit).G2{t)], 
where (Gi(t), G2(t)) possesses joint density 

(5.7) P^G,itlG,it)){s,v) = _^.e-(^^+^^)/(3^/^^), s>0,v>0. 

This result corresponds to (1.7), for i/ = 1/3 and it represents a coun- 
terpart of result (1.6) with the reflecting Brownian motion replaced by the 
product Gi{t) ■ G2(t), with joint distribution given in (5.7). 

In the case u = 2/3^, from (5.1) we have that 

(5.8) ^2/32(2;,*) = —^ / / se-(^'+^')/(=^^)n2/3(x,si;)dsd?; 

2'K\/t Jo Jo 

and this suggests that we interpret U2/32{x,t) as the distribution of the 
process 

(5.9) J2/32{t) = A[Gi{t)-G2{t)]. 
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The process (5.9) is analogous to (5.6) with the role of Brownian motion 
played by the process A. 

Analogously to (5.9), for v = 2/3"^, we get 

'/2/33(t) = A{Gi[Gi{t) ■ G2{t)] ■ G2[Gi{t) ■ G2{t)]} 

which has distribution coinciding with 

U2/33{x,t) 

Q /■+00 /■ + 00 

-3 / ' _(s3_^„3)/(3^3() 



2^T^/t Jo 



2'K\/SV 



se 

+00 c+oo 



«;e-('"'+^')/(=^^)'U2/3(a^, zw) dz dw ) ds dv 



Jo 

as an application of (5.1) and (5.8) shows. 

The results of Theorem 2.1 and 5.1 can be furthermore generalized in 
order to relate the solutions Uy[x,t) with Umu{x,t). 

Theorem 5.2. The solution to equation (2.1), for u <2/m, m>l can 
he represented as 

^{m-l)/2 

Uy[X,t) 



(27r)("^-i)/2Vt 



+ 00 



(5.10) x/ •••/ e-('-r+-+-,';j-i)/ -^^2 

Jo Jo 

X Umuix, WiW2--- Wm.-l) dwi ■ ■ ■ dWm-l- 

Proof. From (1.4), by using the multiplication formula of the Gamma 
function [see Magnus and Oberhettinger (1948)], that is. 



m / \ m J \ m 
= (27r)("^-i)/2^i/2-m^p(^^)^ 

for z = — '^^^^^^ , we get that 
Uy{x,t) 



aWV /;!r(i ^^^+^) ] 



m 



2AW2(27r)(™-i)/2 
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\x\ 



k=0 



\m 2 J \m 2 y 



m 2 



777 

A;!r( l-—u{k + r 



m 



2AW2(27r)(™-i)/2m'"/ 



2v 



■!«m-i,„-l/™„,,-2/m -(m-1)/' 



Wt Wr. ■ ■ -W. 



m—1 



X Umu{x,m '^Wi ■ ■ ■ Wrn~lt) dwi ■ ■■dWm-1- 

By means of the transformation 



we finally get (5.10). □ 

We prove now a general result, valid for any < v <2, which gives another 
representation for the solution u^, = Uu(x,t), alternative to those presented 
in the previous sections. 

Theorem 5.3. The solution to (1.1) with initial condition (1.2) or (1-3) 
has the following form: 

1 r+oo 

u,(x,t) = I e-'"u;^/2-ig-|x|«,''/V{AW^)cosW2) 

^TrXt^' 







(5.11) ''''\-2-^^''\-2))'^'' 

1 f^^ -N.cos(W2)-(Ay)^''''tsin('!^_|^|ysinf— Yuy, 



UTT Jq V 2 V 2 y 

/or < 1/ < 2. 

Proof. By applying the reflection property of the Gamma function we 
rewrite (1.4) as 

(-|x|/(AW2))'= 



'(^'*) 2AW2^fc!r(l 



^^^k\r{i-v/2{k + i)) 
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1 \x\ \' sm{i.7:/2{k + l)) fu 



k=0 ^ ^ ^ 



27rAt^/2 
1 



+00 







(5.12) 

U|^v/2\ 1 e*'""(''+l)/2 _ g-ii'7r(fc+l)/2 



00 



E 



X 
A;=0 

1 



At^/2 / kl 2i 



dw 



27rAt'^/2 



X / e 







2i 



X rg-|x|t«''/2e-/2-/(Ar/2)g»i//27r 



which coincides with the first form of (5.11). The second line can be obtained 
by the change of variable w = {\y)'^/'^t. □ 

Remark 5.2. We can check that, for u = \ (i.e., for the heat equation), 
the first expression in (5.11) reduces to the Gaussian density: 

e '"u;^/^ ^cos ' ,^ dw 



27rAti/2 7o V Ati/2 

(5.13) 

= \w = y\ 



g-x2/(4iA2)_ 



In the last step we used formula 3.896.4, page 514, of Gradshteyn and Ryzhik 
(1994). The same check can be done for the second expression in (5.11). 



An alternative form of (5.11) can be obtained by means of a double inte- 
gration by parts, as follows: 
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Corollary 5.1. The solution to (1.1) with initial condition (1.2) or 
(1.3) can be rewritten as 



Uuix,t) 



(5.14) 



+00 



g-u-g-lxlu;"'' V(Af'^/2) cos(i/7r/2) 



X sm 



X w 



u/2 



sin( — ) ) dw 



for <2. 

Proof. The first integration in (5.11) gives 



Uy{x,t) 



1 



cos 



7ri/|x| sin(i/7r/2) V 2 At''/^ 



VTT \X\W 



v/2 



■ sm 



^ g-«;g-|x|«;'^/2/(Ar/2)cos{i.7r/2) 



+ 



1 



'Kv\x\ sin(z^7r/2) 



+ 00 



g-tUg-|a;|«;'"/V(At'"/^)cos(!/7r/2) 



X COS 



VK \X\W 



u/2 



■ sm 



dw 



+ 



cos(i^7r/2) 
27rsin(i/7r/2)At^/2 Jq 



^-w^u/2-l^-\x\u,-/^/{\t-/-^)cos{u7,/2) 



X cos 



VK \X\W 



v/2 



At'^/2 \2 ) 



sin -— dw 



cot(z^7r/2) 



irvix] 



+ 



TTu\x\ sin(i/7r/2) 



g-iOg-|a;|u;''/V(Ar/2)cos{!/7r/2) 



X cos 



\x\w 



u/2 



V 2 At'^/2 

COs(;/7r/2) |x|«;'-/V(Ar/^)cos(t/7r/2) 



■ sm 



2 



— Uw; 



. , e e 

7rz^|x| sin (z/7r/2) 

UTT \x\w^/'^ 



X sm 



2 At''/2 

cos(z^7r/2) 
'Kv\x\ sin2(z/7r/2) 



■ sm 



g-uig-|a:|u;''/V{At''/2)cos(i/7r/2) 
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X sin 



Z^TT \x\w'^/'^ ( V-k\\ 

7-- sm — aw 



cos2(i/7r/2) 



27rsin2(i/7r/2)AW2 ]q 



g-u,^i./2-lg-|x|«;''/2/(AW2)cos(i/7r/2) 



X sm 



VK \x\w 



vl2 



■ sm 



z^vr 



Therefore, from (5.11) we have that 
cos2(z^7r/2) 



1 + 



sin2(z^7r/2) 
1 



Uy{x,t) 



7rz/|x| sin(z^7r/2) 

"+oo 



g--!«g-|x|io''/2/(At''/2)cos{!/7r/2) 







COS 



i^vr \x\w' 



v/2 



■ sm 



cot 



■Kv\x\ sm(z/7r/2) 



oo 



sm 



VTT \X\W 



v/2 



■ sm 



g -«)g-|x|«;''/2/(At'V2)cos(j/7r/2) 



cos 



sm 



-jcos 



I'TtX 

I sm 



,„ ■ sm — 
Xf^/^ \ 2 



WW 



u/2 



( VTX 

— cot 

V 2 



2 y V 

/ z/vr 



■ sm 



sm 



V2 



cos 



\x\iifl'^ ( vkX 
sin — 



xt^i^ V 2 7 



— COS 



sm 



\x\w 



u/2 



dw. 



dw 



■ sm 



which easily gives (5.14). □ 



Remark 5.3. We can check that, for 1^ = 1, (5.14) reduces again to the 
Gaussian density: 



ui{x,t) 



1 



7r|x| Jo 



e sm 



\x\w 



1/2 



AtV2 



dw 
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x\Jo V2 

-x2/(4tA2) 



\/27r| 



as in (5.13). 

With respect to (5.11), formula (5.14) is more appealing as it allows an 
easier analysis of the limit for — )■ 0: 

. X 1 sin(i/7r/2) f+'^ ,„ , 

(5.15) \m\Uy{x,t) = L \ w'^/^e-'^dw 

1 sin(i/7r/2)^/ 1^^^ 



■KU 



For t —7- +00, (5.15) decreases for all values of ;/ G (0, 2]. 

Moreover in the case 1^ = 1, formula (5.15) gives the maximum value of 
the Brownian density. For v = 2 (5.15) is zero for all t > 0, because in this 
case (1.1) becomes the wave equation and its solution has the form of the 
sum of Dirac's impulse functions travelling in opposite directions. 

By means of the following formula 



f 

Jo 



smqx q 

e ^ ax = arctan -, p> \) 

X p 



[Gradshteyn and Ryzhik (1994), formula 3.941.1, page 523] we can check 
that (5.14) integrates to one, as follows: 



+ 00 o /" + 00 f + OO 



u,{x,t)dx = — / / ie---/V(AW2)cosW2) 

TTZ^ ./o Jo X 



X sin I y^^yj^ ^™ ( ~2~ 



TTU 2 Jq 



■oo 

e-'"dw = l. 



Finally it is interesting to analyze the behavior of the solution (for x 
varying and t fixed), which is substantially different in the two intervals 
< < 1 and 1 < z/ < 2 (see Figure 2 above). We rewrite formula (5.14) as 
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(Al (B) 

Fig. 2. (A) The solution u^{x,t), for <u <1. (B) The solution Uv{x,t), for 1 <u <2. 



follows: for X > 0, 



1 r 

Uy{x,t) = / 



X 



where g{x,w,t) = e'''^^''^^''''/^^ sm{xAsm{vT: /2)) <md A = w"/"^ / Xe/"^ . 
The first derivative of ^^^^'^^ with respect to x is equal to zero if 

g x' 



(5.16) 
where 



+ A sin ^e---^™^('^-/2) cos [xA sin ^) 

= Ae-^^-(-/2)sin(^-xAsin^). 
The solution to (5.16) is 

\gg = \gx + const 

or, otherwise, 

g = xconst. 

By choosing const = 1, we obtain that Uy{x,t) attains its maximum on the 
positive half-line if 

(5.17) xe'^''°''^''^l'^^ =sin{xAsiu'^^. 

For \ <u <2 there exists only one value of x which verifies the condition 
(5.17) and this is in accordance with the behavior of the solutions Uj^ pre- 
sented in Fujita (1990), where the relationship with stable laws is exploited. 
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On the other hand, for < < 1, no positive value satisfies (5.17) and 
therefore the maximum is in the origin. The previous results are confirmed 
by the following theorems. 

We now present the general results concerning the relationship between 
the solution u^{x,t) and the stable densities. We need to analyze the two 
intervals < < 1 and 1 < v <2 separately. 

Theorem 5.4. For < < 1, the solution to 

( d^u _ ^d'^u 
(5.18) {~d¥~ 3;EM,t>0, 

I u{x, 0) = (5(x), 

can he represented as 

1 A^/^i f X^/^t V 



(5.19) 



\ V 1 

where P|(-; f , 1) is the density of a stable distribution of parameters 7 = f 
and T] = 1; see (4-20). 

Proof. From (1.4), by using the reflection formula for the Gamma func- 
tion we have that 

k 



1 / jxl \ 



, AWV kW(l-u(k + l)/2) 

k=0 

(5.20) 

1 ^/ yisin(7ri/(A; + l)/2) /i/(fc + l) 



In view of the series representation of stable functions, which for < a < 1 
reads 



oo 

■ vr 



p„(x;7,l) = -> -1 j X ^ + ' sm 

r=0 



-(7 + a)(r + l) 



[see Feller (1971), formula (6.10), page 583, with some corrections, Lukacs 
(1969) and Zolotarev (1986)], we can obtain the first expression in (5.19). 
The second expression can be derived by applying the self-similarity property 
of the stable random variables. □ 



Finally we consider the case 1 < z^ < 2 and we state the following result: 
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Theorem 5.5. The solution to 
ut{x,0) = 0, 

for l<u<2, can he represented as 



(5.22) 



where P2/u{'\ ~ 1) ^-^ density of a stable distribution of parameters 
7 = ^(z^ — 1) and r] = 1. 

Proof. By following the same steps as in the previous theorem we can 
recognize in (5.20), up to the normalizing constant, the series representation 
of the stable law p2/,y of order a = 2/^ [see (4.21)], so that we get (5.22). 
□ 

Remark 5.4. In view of Theorems 2.3 and 5.5 and by considering the 
property of self-similarity of the stable laws, we can write that 

1 f°° e-^'V(4«'A) I f\w\l \ 

'^ = -Jo ^J^^""'" I AF' ^^2^ - V 

(5.23) 

roc g— z^/(4iuA) / Y 



Pi/^^|u;|,i(2z^- l),A^/''t^ dw. 



V Jo V^TTWX 

Formula (5.23) shows that the solution Ui^, for ^ < < 1, can be inter- 
preted as the distribution of the process B{\Su{t)\), t > 0, where is the 
stable process with density ^Pi/,^{\ ■ |, ^(2z^ — 1), A^/'^t). 

Moreover, as a consequence of Theorems 2.1 and 5.5, the solution of our 
problem (1.1)-(1.2), for i < < 1, can be written in an alternative to the 
form (2.2) also as a stable law evaluated at a Brownian time: 

1 e~''^/(''*) 1 f \x\ I \ 

'^ = ui -1/^2X^P'/^ [X^^ u^^' - V 

Remark 5.5. We check that, for = 1, both expressions (5.19) and 
(5.22) yield the Gaussian density 

(5.24) ui{x,t) = —^e--''^^'^''l 
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We start by considering the last expression in (5.19) 

, , 1 /111 
(5.25) ux(x,t) = --^vxl2\ 



where [from (4.20)], for y > 0, 
(5.26) 

I g-l/(2y(V2Ati/2)2) 

~ V2Ati/2 ■ 

By taking in (5.26) y = we get from (5.25) the Gaussian density (5.24). 
Formula (5.22) immediately supplies (5.24) for v = I. 

Acknowledgments. The authors thank one anonymous referee for bring- 
ing their attention to some relevant papers on fractional equations. Thanks 
are also due for his accurate check of the text and of the calculations. 
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